Quantum simulator for the 0(3) nonlinear sigma model 
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We propose a design for the construction of a laboratory system based on present-day technology 
which reproduces and thereby simulates the quantum dynamics of the 0(3) nonlinear sigma model. 
Apart from its relevance in condensed-matter theory, this strongly interacting quantum field theory 
serves as an important toy model for quantum chromo-dynamics (QCD) since it reproduces many 
crucial properties of QCD. 



PACS numbers: 03.67.-a, 03.67.Lx, ll.lO.Kk, 68.65.-k. 
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Introduction In many areas of physics, progress has 
been thwarted by our lack of understanding strongly in- 
teracting quantum systems with many degrees of free- 
dom such as quantum field theories. Beyond perturba- 
tion theory with respect to some parameter or semiclassi- 
cal models/methods, there are not many analytical tools 
available for the treatment of these systems. Numeri- 
cal methods are hampered by the exponentially increas- 
ing amount of resources required for the simulation of 
quantum systems with many degrees of freedom in gen- 
eral. However, this obstacle applies to classical comput- 
ers only - quantum computers will be able to simulate 
other quantum systems with polynomial effort I*]. But 
as long as universal quantum computers of sufhcicnt size 
(e.g., number of QuBits) are not available, one has to 
search for alternatives. One possibility is to design a 
special quantum system in the laboratory which repro- 
duces the Hamiltonian of a particular quantum field the- 
ory of interest. This designed quantum system can then 
be regarded as a special quantum computer (instead of a 
universal one) which just performs the desired quantum 
simulation. In this Letter, we propose such a quantum 
simulator for the example of the 0(3) nonlinear sigma 
model and demonstrate that it can be constructed using 
present-day technology. 

The Model The l+l dimensional 0{N) cr-model % 
H H H 0, 0, H is described by the 0{N) and Poincarc 
invariant action 
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with the internal vector cr = (cti, (72, ... , cfn) G re- 
flecting the 0(A^)-symmetry. So far, this theory describes 
independent free fields, but the constraint 



'N 



N 

7^2 



(2) 



introduces an interaction corresponding to the coupling 
5 > 0. For vanishing coupling g [ the curvature of 
the constraint sphere (cr^ = N/g^) vanishes and we re- 
produce (locally) an effectively free theory. For finite 
coupling (7 > 0, we obtain a non-trivially interacting the- 
ory as long as iV > 3, i.e., in the non-Abelian case. The 



classical ground state cr = const is 0(iV)-degenerate, but 
quantum interaction lifts that degeneracy and gives the 
classical Goldstone modes a mass gap, see, e.g., 

Properties Apart from its relevance in condensed- 
matter theory (partly due to its relation to spin-systems 
such as anti-ferromagnets, see, e.g., 0), it can be shown 
that the 1-1-1 dimensional 0{N) ct- model reproduces 
many crucial properties of quantum chromo-dynamics 
(QCD) and hence serves as an important toy model, 
see, e.g., 's^. The cr-model is renormalizable (in l-fl 
dimensions, cf. ^]) and its running coupling g{p^) gen- 
erates asymptotic freedom g^{p^ 3> A'^) cx l/ln(p^/A^), 

jTTL 

cf. 3]. In analogy to QCD, the classical scale invariance 
— !■ Q.x^ is broken dynamically corresponding to the 
dimensional transmutation 5 — )■ A. Furthermore, the a- 
modcl generates non- vanishing vacuum condensates such 
as (£) 7^ in the operator product expansion and repro- 
duces the trace anomaly {T^) 7^ (see, e.g., 0). It also 
serves as a toy model for the study of the low-energy 
theorems and sum rules (see, e.g., |3). For iV = 3, 
the cr-model exhibits instantons (mapping of §2 onto K^, 
cf. [IQ). Furthermore, the model is exactly solvable in 
the large- A'^ limit, where it corresponds to massive free 
fields with sub-leading (in 1/A^) interaction terms, see, 
e.g., Q. Its complex version, the CP{N — l)-model, re- 
produces confinement (though this is not such a strik- 
ing feature in 1-1-1 dimensions) and the J7(l)-problcm of 
QCD (see, e.g., i). 

The Analogue In order to reproduce the quantum dy- 
namics of the the H-l dimensional 0{N) cr-model accord- 
ing to Eqs. 1^ and (O, let us consider a large number of 
perfectly insulating thin hollow spheres with the radius 
p lined up at equal distances Ax with single electrons 
being captured by the polarizability (inducing a finite 
extraction energy) on each of the hollow spheres. These 
insulating spheres are surrounded by an arrangement of 
superconducting spheres (radius a) and wires (radius S) 
as depicted in Fig. ^ which generate controlled interac- 
tions of the confined electrons via their image charges. 
The involved length scales including the typical wave- 
length of the excitations A, the distance of elements (lat- 
tice spacing) Ax, the distance between the insulating and 
the conducting spheres 7, the radii of the insulating and 
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conducting spheres p and a and wires 5 (cf. Fig. ^ are 
supposed to obey the foUowing hierarchy 

(3) 



A > Ax > 7 > p, a 
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FIG. 1: Sketch of the proposed analogue quantum simulator. 
The solid lines and spheres denote (super) conductors and 
the hollow spheres are insulators containing single electrons. 
Shown are just three elements of a long chain (top) and a 
close-up view (bottom) with the involved length scales. 

The total Lagrangian for the system of electrons reads 



vifith m being the mass of the electrons and V{ri+i,ri) 
their interaction potential, where only nearest neighbors 
arc taken into account in view of the assumptions 0. In 
this limit, the interaction potential induced by the images 
of the electron charges e simplifies to 



F(ri+i,r,) 
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where the first addend in the denominator on the right- 
hand side is due to the capacitance of the conducting 
spheres Aneoct and the second one due to the capacitance 
of the long wires 2TTeoAx/\n{Ax/S). Comparing the re- 
sulting Lagrangian in Eqs. I0J and © with the one in 
Eq. we can read off the effective propagation speed 



Coff = Co 4 



Aneomcl 4a + Ax/ \ii{Ax/S) 
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Since the first term under the root represents the classical 
electron radius (of order 10^^^ m), the effective propa- 
gation speed Ccff is much smaller than the speed of light 
in vacuum cq ^ CqB for realistic parameters (see below), 
i.e., we obtain a large slow-down. Furthermore, we may 
identify the effective coupling for iV = 3 
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where the first term under the root is the classical elec- 
tron radius over the square of the fine structure constant. 
The value of the effective coupling can be tuned by vary- 
ing the ratio 7/p 3> 1 and may well be of order one (see 
parameters below). Strictly speaking, the above equation 
determines the value of the running coupling (7eff(p^) at 
a length scale corresponding to the lattice spacing Ax 
(lattice renormalization scheme). In complete analogy to 
Aqcd, the coupUng gl^(p^ > A^) oc \/\n{p^/Kl) deter- 
mines the induced scale of dynamical symmetry break- 
down Act of the (T- model (dimensional transmutation). 
This important quantity sets all other length scales such 
as the mass gap (see, e.g., 0) and must satisfy the con- 
dition (PJ for consistency, i.e., K^Ax <C 1. Finally, iden- 
tifying (again for = 3) 



a{x = lAx) = 
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(8) 



the continuum limit (^ • Aa; — s- ^ dx for A 3> Ax) of 
Eq. Q generates the Lagrangian Q of the 0(3) non- 
linear sigma model with the constraint |(5J) being imple- 
mented by rf = 

Disturbances Of course, for a realistic proposal, it 
is essential to estimate the impact of the contributions 
which have been omitted so far. The additional kinetic 
terms due to inductances L of the wires are negligible 



(4) LI <C mf provided that 
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holds, i.e., for a sufficiently large slow-down (as one would 
expect). For the same reason, the influence of the zero- 
point fluctuations of the electromagnetic field (induc- 
tance of free space) is negligible. 

In contrast to sequential quantum algorithms, where 
errors may accumulate over many operations, the quan- 
tum simulation under consideration is basically a ground 
state problem and hence more similar to adiabatic quan- 
tum computing |lOj . In this case, decoherence can be 
neglected as long as the interaction energies of the dis- 
turbances are much smaller than the energy gap between 
the ground state and the first excited state |lO|. For 
the nonlinear (T-model, this gap is determined by the in- 
duced scale Act (in analogy to QCD). Therefore, the ener- 
gies of all perturbations (e.g., impurities in the material) 
must be much smaller than the gap of order hccs^a- In 
particular, in order to see quantum behavior (where the 
Heisenberg uncertainty relation becomes important), the 
temperature must be small enough 



fceT < hccsAo 



(10) 



Another issue concerns the spins of the electrons, which 
have been omitted so far. Fortunately, we may fix the 
electron spins by a small external magnetic field (see the 
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next paragraph) and the various spin-spin and especially 
spin-orbit coupling terms are negligible (in comparison 
to fiCeffAcr) for the parameters provided below. 

Phase Diagram Before investigating the impact of an 
external magnetic field, let us turn to the phase dia- 
gram of the nonlinear a-model in terms of the temper- 
ature T and the chemical potential /i. For low tem- 
peratures k^T <C hccs^a- and small chemical potentials 
fi ^ SceffAo-, we basically get the usual vacuum state. 
Note that the introduction of a chemical potential ne- 
cessitates the definition of a particle number (which is a 
nontrivial issue in interacting theories). In the ct- model, 
this can be achieved by means of the Noether current cor- 
responding to the global 0(3) invariance = <t x d^cr 
and the associated global charge along some internal axis 
n with = 1 



Q — n ■ I dx cr X (T . 

Ccff 



(11) 



For the laboratory system, the Noether charge Q is just 
the total (orbital) angular momentum in units of h. Note 
that still many charges Q :s> I are required to generate 
one magnetic flux quantum (due to cq ^ Ccff). 

In terms of the chemical potential defined with re- 
spect to this (dimensionless) Noether charge, the grand- 
canonical Hamiltonian H^c reads 



(12) 



Translating this expression back to our laboratory system 
in Eq. Q , we observe that the chemical potential exactly 
corresponds to an external magnetic field B inducing the 
additional term r ■ A = r ■ {r x B)/2> ~ B ■ {r x r)/3 
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(The second-order term e^A^jm is three orders of mag- 
nitude smaller for the parameters given below and can 
be neglected.) When the effective chemical potential ^cff 
exceeds the energy gap of order ficesAo-, the structure 
of the ground state changes and the above Noether cur- 
rent 2i, acquires a non- vanishing expectation value (quan- 
tum phase transition, see, e.g., 0). At the critical field 
Scrit = C(™coftAcr/e) where this quantum phase tran- 
sition occurs, the energy of the electron spins is of the 
same order as the gap fi^ ■ B = 0{hcceA^) and thus 
much bigger than the temperature. Hence one can fix 
the electron spins with much smaller external magnetic 
fields B <C i?ciit without disturbing the vacuum state too 
much. On the other hand, it is also possible to explore 
the full phase diagram (e.g., cross the quantum phase 
transition, monitored by a SQUID) by increasing the ex- 
ternal magnetic field - which is completely equivalent to 
changing the chemical potential (and hence the number 
of particles). For the set of parameters discussed below, 
the critical field Bcrit is of order milli-Tesla. 



Experimental Parameters The aforementioned con- 
straints, in particular Eqs. (O and H10(l . provide the 
frame of a window of opportunity for the experimental 
realization of the proposed quantum simulator - which is 
(fortunately) open to present-day technology. Let us first 
explore the limit set by the ultra- low temperatures. For 
solid bodies of reasonable size, one can reach temper- 
atures of order 10 /iK by electron gas cooling via spin 
relaxation. If we choose our parameters according to 
S = 100 nm, p = 400 nm, a = 500 nm, 7 = 2.5 /im, and 
Ax — 12.5 fim, we obtain gcs = 0{1), CcS ~ lO^m/s, 
A~^ « 125 ^m, and hccs^a corresponds to 600 /iK, 
which satisfies all of the above assumptions reasonably 
well. Alternatively, we may start from the present state 
of nanotechnology which facilitates the production of 
nanowircs with a radius of order nanometer. If we explore 
this limit and choose S — 1 nm, p = 12 nm, a — 5 nm, 
7 = 25 nm, and Aa; — 125 nm, we obtain a similar value 
for g^s and Ccs ~ lO^m/s, but now hccsAa- corresponds 
to a temperature of order Kelvin. The range between 
/iK and fractions of a Kelvin as well as between nanome- 
ters and micrometers provides a two or three orders of 
magnitude wide window of opportunity and the optimum 
experimental parameters are probably somewhere in the 
middle. 

The thin superconducting wires can be switched on 
and off by local variations of the temperature (below and 
above the critical value). If the interaction V{ri^i,r.i) 
is switched off, the energy spectrum of the electrons is 
determined by the usual spherical harmonics 



h'^ £{£+!) 
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The energy gap between the s-state ^ = and the p-state 
£ = 1, i.e., without interaction V{ri+i,ri), is one order 
of magnitude larger than with interaction hces^a-- Con- 
sistently, the interaction potential V{ri+i,ri) between 
the electrons on different spheres is of the same order of 
magnitude as the gap between the s-state £ = and the 
p-state ^ = 1 on a single sphere leading to strong entan- 
glement of the ground state. If we want to switch on the 
interaction V{ri+i,ri) adiabatically (e.g., via changing 
the temperature of the wires) satisfying the condition for 
the adiabatic theorem | {ipo\ dH{t)/dt \tl^i) \/{AEoif < 1 
in order to stay in the ground state |V'o), the typical adi- 
abatic switching time should be longer than a few pi- 
coseconds. Finally, for the parameters discussed above, 
the various spin-spin and spin-orbit coupling energies are 
at least two orders of magnitude smaller than hccs^a- 

Summary As we have demonstrated above, it is pos- 
sible to construct a quantum simulator for the 0(3) 
nonlinear cr-model with present-day technology. Such a 
restricted quantum computer would allow the compar- 
ison, for a controllable scenario, between perturbative 
and non-perturbative analytical methods (renormaliza- 
tion ffow 0, H, , instantons [3| , operator product ex- 
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pansion and vacuum condensates, low-energy theorems 
and sum rules [5|, the S-matrix [6| etc.) as well as nu- 
merical results |3 on the one hand with real quantum 
simulations on the other hand. In contrast to most of the 
numerical simulations, for example, the proposed quan- 
tum simulator works in real (laboratory) time, i.e., it is 
not necessary to perform a Wick rotation to Euclidean 
time. This advantage facilitates the study of the evolu- 
tion of excitations, for example collisions (S-matrix etc.). 

Furthermore, the proposed set-up allows a direct ac- 
cess to the quantum state and hence an investigation of 
the strong entanglement (e.g., in the ground state or near 
the quantum phase transition). This could be done via 
state-selective radio/micro- wave spectroscopy of transi- 
tions from the levels in Eq. H14(l to some higher-lying 
empty and isolated internal level (of the semi-conductor) 
with a sharp energy, for example (fluorescence measure- 
ment). Generating the radio/micro- waves via a circuit 
(wave-guide) facilitates the position control of the mea- 
surement (vicinity of the inductance loop). Furthermore, 
one may also switch off the wires (e.g., by locally increas- 
ing the temperature) before the measurement. 

It is also possible to create particles (and their anti- 
particles), which can be used to study the S'-matrix, for 
example, via the illumination with (left and right) cir- 
cular polarized radio/micro- wave radiation, cf. Eq. 
and the subsequent remarks. Another interesting point is 
the robustness or fragility of non-perturbative properties 
(such as the instanton density) with respect to a small 
coupling to external degrees of freedom. 

Apart from above points of interest, the construction 
of such a restricted quantum computer, which is espe- 
cially dedicated to the simulation of the 0(3) nonlin- 
ear cr-model, would be an interesting feasibility study for 
more general quantum simulators for a comparably well 
understood (yet nontrivial) system. Finally, experience 
shows that the availability of a new tool (such as the 
proposed quantum simulator) yielding new tests/results 
usually leads us to a new level of understanding in physics 
with possibly unexpected outcomes. 

Outlook After having handled and understood the 1+1 
dimensional situation, the extension to 2+1 dimensions 
should not be very problematic. The 2+1 dimensional 
0(3) nonlinear cr-model loses some of the properties dis- 
cussed above, but also acquires novel features, such as 
skyrmions which are described by the topological cur- 
rent jP = e'^'^P (T ■ {d^cr X d^cr). The inclusion of an 
explicit 0(3)-symmetry breaking term n ■ cr should be 
easy in 1+1 and 2+1 dimensions. Note that we did not 
incorporate a topological (Chern-Simons type) ^-term 
Ce = e^'^ <j ■ {dfj^cr x di^cr) in our 1+1 dimensional sce- 
nario (in analogy to the ^^-ternl G*^G''"' in QCD), whose 
implementation is less straight-forward. Further interest- 
ing topics are the behavior of strongly interacting quan- 



tum field theories (such as QCD and the cr-model) during 
the cosmic expansion and the (long-range) entanglement 
of QCD vacuum state (which might be used as a tool for 
diagnosis and a resource). 
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